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Abstract. The antisymmetrized molecular dynamics (AMD) simulations 
suggest that the isospin composition of fragments produced dynamically in 
multifragmentation reactions is basically governed by the symmetry energy of 
low-density uniform nuclear matter rather than the symmetry energy for the 
ground-state finite nuclei. After the statistical secondary decay of the excited 
fragments, the symmetry energy effect still remains in the fragment isospin 
composition, though the effect in the isoscaling parameter seems a very delicate 
problem. 
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1. Introduction 

In medium energy heavy-ion collisions, many fragments are formed almost simulta- 
neously in an expanding and excited nuclear system, which gives us an opportunity 
to investigate the nuclear matter properties at various temperatures and densities, 
such as the nuclear equation of state and the liquid-gas phase transition. However, 
fragments are formed in dynamically evolving system in most cases, and therefore 
dynamical model calculations are necessary. 

Antisymmetrized molecular dynamics (AMD) model |2 El EI E] respects 
several quantum features in fragment formation reactions. It uses a fully antisym- 
metrized many-body wave function of Gaussian wave packets, which can describe 
the ground state properties of nuclei reasonably well. The single-particle evolution 
in the mean field potential is described by the the motion of the wave packet cen- 
troids and the stochastic quantum branching process which respects the change of 
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the shape of the phase-space distribution, keeping the advantage of the molecular 
dynamics that the fragments are formed with the wave packets localized in phase 
space. The two-nucleon collision effect is also treated as a stochastic process. 

The AMD simulations for heavy-ion collisions are useful not only to explain 
the experimental data but also to know what kind of information is reflected in the 
fragment formation. In Refs. [El E]> we have analyzed the fragment yields in the 
AMD simulations for multifragmentation reactions of the central collisions of Ca 
isotopes at 35 MeV/nucleon, in order to see how the fragment isospin composition 
is related to the symmetry energy term of the effective interaction adopted in the 
calculation. Isoscaling has been observed in the fragment yield ratios from two 
reaction systems with different proton-to-neutron ratios in the AMD results as well 
as in the experimental data [ and in the predictions by various statistical models [ 
ITT)] and other dynamical models [E] El The symmetry energy has turned out 
to be reflected in the fragment isospin composition with almost no surface effect, 
which means that the property of low-density uniform nuclear matter is responsible 
for the fragment isospin composition rather than the symmetry energy for finite 
nuclei. The results for the primary fragments are discussed in Sec. 2. 

The fragments are recognized in the AMD simulations at a finite time t ~ 300 
fm/c. These primary fragments are excited with the typical excitation energy of 
about 3 MeV/nucleon and they will decay by emitting particles with a long time 
scale before they are detected in experiments. In order to compare the calculated 
results with the experimental data, the secondary decay of primary fragments should 
be considered by employing a statistical decay code. The effect of the secondary 
decay on the symmetry energy effects is studied in Sec. 3. Some features observed 
in the primary fragments remain after the secondary decay, which can be utilized to 
get the information of the symmetry energy from the experimental data in principle. 
However, the symmetry energy effect in the isoscaling parameter is produced by a 
very delicate cancellation of the effect in the width of the isotope distribution and 
that in the mean isospin asymmetry of fragments. 

2. Symmetry Energy Effects in Primary Fragments 

In this section, we discuss how the density-dependent symmetry energy is reflected 
in the fragment isospin composition, by using the result of AMD simulations of 
Refs. [El E| for 40 Ca + 40 Ca, 48 Ca + 48 Ca, and 60 Ca + 60 Ca. We simulate collisions 
by boosting two nuclei whose centers are separated by 9 fm and calculating the 
dynamical evolution of each collision until t = 300 fm/c. The primary fragments 
are recognized at t = 300 fm/c by the condition that the two nucleons belong to 
the same fragment if the spatial distance between them is less than 5 fm. The 
calculations are done with the Gogny force [E] and the Gogny-AS force [El which 
are different in the density dependence of the symmetry energy for nuclear matter 
as shown in Fig. ^ 

We have seen [Ej that the fragment yield ratios Y2(N,Z)/Yi(N,Z) between 
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Fig. 1. Density dependence of the symmetry energy of nuclear matter for the 
Gogny force and for the Gogny-AS force. 

different reaction systems satisfy isoscaling 

Y 2 (N,Z)/Y 1 (N,Z) oc e aN+0z . (1) 

Isoscaling is equivalent to the fact that the fragment yields of different reaction 
systems i are given by 

Y i (N,Z)=exp[-K(N,Z)+a i N + i Z + <y i ], (2) 

where fli and 7$ are constants that depend on the reaction system i, while 
K(N, Z) is a function that is independent of the reaction system. By combining the 
fragment yields from the three reaction systems, we can get the function K(N, Z) 
for a wide region of [N, Z) , even though the number of generated events is not very 
large. 

The obtained K(N,Z) behaves very smoothly as a function of N and Z [Ej. 
The shell and paring effects are weak in K(N, Z) compared to the effects in the 
ground state binding energies. We find [[7] that K(N, Z) can be fitted very well by 
the functional form 

K(N, Z) = £(Z)N + rj(Z) + C{Z) ^ (3) 

for each Z . The quantity C(Z) is related to the width of the isotope distribution 
for the given and it is expected to be sensitive to the symmetry energy. 

The top panel of Fig.|21shows the obtained parameter Q{Z) as a function of the 
fragment proton number Z for the Gogny force and the Gogny-AS force. We find 
that the Z-dependence of ((Z) is very weak for Z > 5, which is not consistent with 
the idea of the usual equilibrium of primary fragments. In fact, if the fragments 
were in thermal and chemical equilibrium, C(Z) would be related to the symmetry 
energy for the fragment nuclei as (,(Z) = (c v + c s A _1 / 3 )/T, where the volume and 
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surface symmetry energy coefficients satisfy c v ~ — c s for the ground state nuclei. 
This size dependence of the symmetry energy is not compatible with the very weak 
Z-dependence of Q{Z) for Z > 5 [0. 

It may be rather surprising that we can find almost no surface effect even 
though we are looking at relatively small fragments. There can be several possible 
explanations [u\. One of them is associated with the fact that fragments are not 
isolated when they are formed. When the density fluctuation is developing from a 
uniform low density matter, the fragments are still interacting with attractive force 
through their surfaces. Therefore, the surface free energies can be expected to be 
smaller for these fragments than for the isolated fragments. Independent of the 
physical origin for the weakening of the surface symmetry free energy, it suggests 
that the volume quantity, which is the same as that in the infinite nuclear matter, 
can be directly obtained by the analysis of the fragmentation results even though 
the produced fragments are not very large. 

If we adopt the interpretation that £ is related to the symmetry energy C sym (p) 
for the uniform matter at a certain density p by 

C=%^, (4) 

we can utilize the symmetry energy effect observed in the top panel of Fig. [21 to 
derive p and T [ Ej • In order to explain the ratio of £ for the Gogny force and 
the Gogny-AS force which are different in the density dependence of the symmetry 
energy as shown in Fig. ^ the density should be p ~ \p§- Furthermore, from the 
absolute values of £ and C sym (ipo), the temperature should be T ~ 3.4 MeV. This 
condition of the density and the temperature is reasonable as the condition for the 
fragment formation, which is necessary for the justification of Eq. (JIJ. 

If the fragment yields from different reaction systems satisfy isoscaling, the 
isoscaling parameter a is related to Q(Z) by [0 

a = ±aZ) x HZ/Af, (5) 

where 

A(Z/A) 2 = (Z/A^Z)) 2 - (Z/A 2 (Z)) 2 (6) 

represents the difference of the mean isospin asymmetry of fragments for each given 
Z between the two reaction systems. We will choose the 40 Ca + 40 Ca system as the 
system 1 in this paper. Equation J5j can be derived from Eq. (0) by employing Eq. 
©, and thus it is not relevant how ((Z) is related to the symmetry energy. It should 
be noted that a and A(Z/A) 2 depend on the reaction systems while ((Z) does not, 
and that (,(Z) and A(Z/A) 2 are functions of Z while a should be independent of 
Z. 

The middle panel of Fig.[21shows A(Z/A) 2 as a function of Z when the 60 Ca + 
60 Ca (and 48 Ca + 48 Ca) system is chosen as system 2. The results with the two 
different symmetry energy terms are shown. The symmetry energy effect is clearly 
seen in A(Z/A) 2 as well as in ((Z), though the effect is in the opposite direction. 
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Fig. 2. The three quantities ((Z) Fig. 3. The same as Fig. |3 but for the 

(top), A(Z/A) 2 (middle) and a (bot- final fragments, 
torn) for the primary fragments. The 
results with the Gogny and Gogny-AS 
forces are shown by filled diamonds and 
the open squares, respectively. 

The bottom panel shows 4^(Z) x A(Z/A) 2 which is the product of the top and 
middle panels and should coincide with the isoscaling parameter a. The product is 
actually almost independent of Z, which is consistent with Eq. Ipjl. The symmetry 
energy effect in a is similar to the effect in ((Z), but it has been produced by the 
cancellation of the effect in ((Z) and the effect in A(Z/A) 2 , with the former being 
larger than the latter. 



3. Symmetry Energy Effects in Final Fragments 

The statistical decay of primary fragments is calculated by using the code [ O] 
based on the sequential binary decay model by Piihlhofer [EI- The code employed 
in the present work also takes account of the emission of composite particles not 
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only in their ground states but also in their excited states with the excitation energy 
E* < 40 MeV. The experimental information is incorporated for known levels of 
A < 28 nuclei, while the Fermi-gas level density is assumed otherwise with the level 
density parameter a — A/(8 MeV). 

For the final fragments after the statistical decay, the quantities C,(Z) and 
A(Z/A) 2 can be defined in the same way as for the primary fragments. Equa- 
tion 10 holds for the final fragments as well if isoscaling is good and Eq. @ is 
approximately valid for the final fragments. However, C,{Z) for the final fragments 
is not necessarily related to the symmetry energy directly in a simple way. 

The top and middle panels of Fig. [3] show £(Z) and A(Z/A) 2 , respectively, 
calculated for the final fragments. The absolute values of (,{Z) and A(Z/A) 2 change 
from those in Fig. for the primary fragments. Nevertheless, the effect of the 
symmetry energy term in (,{Z) and A(Z/A) 2 is similar to what we have observed 
for the primary fragments. This means that the symmetry energy effect in the 
dynamical stage of the reaction is observable in principle even after the secondary 
decay of fragments. 

However, in the isoscaling parameter a which is shown as 4£(Z) x A(Z/A) 2 in 
the bottom panel of Fig. the symmetry energy effect is not clearly seen. This is 
because the symmetry energy effect in C(Z) is largely canceled by that in A(Z/A) 2 
to give this small effect in a. We expect that the situation of such cancellation may 
be sensitive to the details of the statistical decay code, and therefore it is not easy 
to predict the symmetry energy effect in a for the final fragments at the current 
stage. 

4. Summary 

The isospin composition of fragments produced by the AMD simulations is analyzed 
to see how the symmetry energy is reflected. The studied quantities Q{Z), A(Z/A) 2 
and a for the primary fragments are sensitive to the symmetry energy term in the 
effective interaction, and the calculated results suggests that the symmetry energy 
of uniform nuclear matter at a reduced density p ~ is directly responsible for 
the isospin composition of the primary fragments even though the fragments are of 
finite size. A statistical decay calculation shows that the symmetry energy effects 
in (,{Z) and A(Z/A) 2 remain after the secondary decay of the primary fragments. 
However, the symmetry energy effect in a for the final fragments requires a careful 
study because it is given by a delicate cancellation of the effect in C,(Z) and the 
effect in A(Z/A) 2 . 

Acknowledgments 

The main part of this talk was based on the collaboration with M. B. Tsang, W. 
G. Lynch, P. Danielewicz and W. A. Friedman, which was supported by Japan 
Society for the Promotions of Science and the US National Science Foundation 



Multifragmentation and Symmetry Energy Studied with AMD 



7 



under the U.S. -Japan Cooperative Science Program (INT-0124186), by the High 
Energy Accelerator Research Organization (KEK) as the Supercomputer Project, 
and by grants from the US National Science Foundation, PHY-0245009, PHY- 
0070161 and PHY-01-10253, and a Grant-in-Aid for Scientific Research from the 
Japan Ministry of Education, Science and Culture. The work was also partially 
supported by RIKEN as a nuclear theory project. 

References 

1. A. Ono and H. Horiuchi, Porg. Part. Nucl. Phys. 53, 501 (2004). 

2. A. Ono, H. Horiuchi, Toshiki Maruyama and A. Ohnishi, Phys. Rev. Lett. 
68, 2898 (1992); A. Ono, H. Horiuchi, Toshiki Maruyama and A. Ohnishi, 
Prog. Theor. Phys. 87, 1185 (1992). 

3. A. Ono and H. Horiuchi, Phys. Rev. C 53, 2958 (1996). 

4. A. Ono, Phys. Rev. C 59, 853 (1999). 

5. A. Ono, S. Hudan, A. Chbihi, J. D. Frankland, Phys. Rev. C 66, 014603 
(2002). 

6. A. Ono, P. Daniclcwicz, W.A. Friedman, W.G. Lynch and M.B. Tsang, Phys. 
Rev. C 68 051601 (R) (2003). 

7. A. Ono, P. Daniclcwicz, W.A. Friedman, W.G. Lynch and M.B. Tsang, Phys. 
Rev. C 70 041604(R) (2004). 

8. J. Decharge and D. Gogny, Phys. Rev. C21, 1568 (1980). 

9. H. S. Xu, M. B. Tsang, T. X. Liu, X. D. Liu, W. G. Lynch, W. P. Tan, A. 
Vander Molen, G. Verde, A. Wagner, H. F. Xi, C. K. Gelbke, L. Beaulieu, B. 
Davin, Y. Larochelle, T. Lefort, R. T. de Souza, R. Yanez, V. E. Viola, R. J. 
Charity and L. G. Sobotka, Phys. Rev. Lett. 85, 716 (2000). 

10. M. B. Tsang, C. K. Gelbke, X. D. Liu, W. G. Lynch, W. P. Pan, G. Verde, H. 
S. Xu, W. A. Friedman, R. Donangelo, S. R. Souza, C. B. Das, S. Das Gupta, 
and D. Zhabinsky, Phys. Rev. C 64, 054615 (2002). 

11. TX. Liu, X.D. Liu, M.J. van Goethem, W.G. Lynch, R. Shomin, W.P. Tan, 
M.B. Tsang, G. Verde, A. Wagner, H.F. Xi, H.S. Xu, M. Colonna, M. Di 
Toro, M. Zielinska-Pfabe, H.H. Wolter, L. Beaulieu, B. Davin, Y. Larochelle, 
T. Lefort, R.T. de Souza, R. Yanez, V.E. Viola, R.J. Charity, L.G. Sobotka, 
Phys. Rev. C 69, 014603 (2004). 

12. M. Colonna, F. Matera, Phys. Rev. C 71, 064605 ( 2005). 

13. C. O. Dorso, C. R. Escudero, M. Ison, J. A. Lopez, |nucl-t h/0504036| 

14. Toshiki Maruyama, A. Ono, A. Ohnishi and H. Horiuchi, Prog. Theor. Phys. 
87, 1367 (1992). 

15. F. Puhlhofer, Nucl. Phys. A280, 267 (1977). 



